Quantum transport of two-dimensional Dirac fermions in SrMnBi2 
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We report two-dimensional quantum transport in SrMnBi2 single crystals. The linear energy dis- 
persion leads to the unusual nonsaturated linear magnetoresistance since all Dirac fermions occupy 
the lowest Landau level in the quantum limit. The transverse magnetoresistance exhibits a crossover 
at a critical field B* from semiclassical weak-field B^ dependence to the high-field linear-field de- 
pendence. With increase in the temperature, the critical field B* increases and the temperature 
dependence of B* satisfies quadratic behavior which is attributed to the Landau level splitting of the 
linear energy dispersion. The effective magnetoresistant mobihty fiMR 3400 cm^/Vs is derived. 
Angular dependent magnetoresistance and quantum oscillations suggest dominant two-dimensional 
(2D) Fermi surfaces. Our results illustrate the dominant 2D Dirac fermion states in SrMnBi2 and im- 
ply that bulk crystals with Bi square nets can be used to study low dimensional electronic transport 
commonly found in 2D materials like graphene. 
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Dirac fermions have raised great interest in condensed 
matter physics, as seen on the example of materials such 
as graphenei and topological insulators (TIs)<^ The hn- 
ear dispersion between momentum and energy of Dirac 
fermions brings forth some spectacular properties, such 
as zero effective mass and large transport mobilityji— 
In addition to the surface/interface states in TIs and 
graphene, Dirac states in bulk materials were discussed 
in organic conductors^ and iron-based superconductors 
such as BaFe2As2<^'^ Recently, highly anisotropic Dirac 
states were observed in SrMnBi2 j^*^ where linear energy 
dispersion originates from the crossing of two Bi Gpx.y 
bands in the double-sized Bi square nets. SrMnBi2 
has a crystal structure similar to that of the su- 
perconducting Fe pnictides and is a bad metal^i^ 
The Fermi velocity along T — M symmetry line is Vp « 
1.51 X 10^ m/s, whereas the Fermi velocity in the orthog- 
onal direction experiences nearly one order of magnitude 
decreaseiii^ 

One of the interesting properties of Dirac materials is 
the quantum transport phcnomcnai^iiS Unlike the con- 
ventional electron gas with parabolic energy dispersion, 
where Landau levels (LLs) are equidistant the distance 
between the lowest and 1'** LLs of Dirac fermions in mag- 
netic field is very large and the quantum limit where all 
of the carriers occupy only the lowest LL is easily realized 
under moderate fieldsJ^i^ Consequently some quantum 
transport phenomena such as quantum Hall effect and 
large linear magnetoresistance (MR) could be observed 
by conventional experimental methods in Dirac fermion 
systemjii"— 

Here we show two-dimensional (2D) quantum trans- 
port in bulk SrMnBi2 single crystals. The linear energy 
dispersion leads to the unusual nonsaturated linear MR 
since all Dirac fermions occupy the lowest LL in the quan- 
tum limit. The transverse MR exhibits a crossover at a 
critical field B* from semiclassical weak-field MR ^ to 



the high-field MR ~ B dependence. The critical field B* 
increases with the increase in temperature and its tem- 
perature dependence satisfies quadratic behavior which 
is attributed to the Landau level splitting of the linear 
energy dispersion. Angular dependent MR and oscilla- 
tion indicates the quasi-2D Fermi surfaces (FS). We de- 
rive the effective magnetoresistant mobility hmr ^ 3400 
cm^/Vs. Our results illustrate the dominant 2D Dirac 
fermion states in SrMnBi2. 

Single crystals of SrMnBi2 were grown using a self- 
flux methodfiS Stoichiometric mixtures of Sr (99.99%), 
Mn (99.9%) and excess Bi (99.99%) with ratio 
Sr:Mn:Bi=l:l:9 were sealed in a quartz tube, heated to 
1050 °C and cooled to 450 °C where the crystals were 
decanted. X-ray diffraction (XRD) data were taken with 
Cu Ka (A = 0.15418 nm) radiation of Rigaku Miniflex 
powder diffractometer. Transport measurements were 
conducted in a Quantum Design PPMS-9 with conven- 
tional four-wire method. The crystal was cleaved to a 
rectangular shape with dimension 4x1 mm^ in the ab- 
planc and 0.2 mm thickness along the c-axis. For in-plane 
resistivity pab(T), the current path was in the a&-plane, 
whereas magnetic field was perpendicular to the current 
and parallel to the c-axis except in the rotator exper- 
iments. The c-axis resistivity Pc{T) was measured by 
attaching current and voltage wires to the opposite sides 
of the plate-like crystals High field MR oscillation were 
performed at National High Magnetic Field Laboratory 
in the same configuration to the in-plane MR. 

All powder and single crystal XRD refiections can be 
indexed in the I4/mmm space group by RIETICA soft- 
ware (Fig. l(a))i^ The determined lattice parameters 
are a = b = 0.4561(8) nm and c = 2.309(6) nni in agree- 
ment with the published data^^ The in-plane resistivity 
Pab{T) shown in Fig. 1(b) exhibits a metallic behavior. 
An external magnetic field enhances the resistivity. As 
the temperature is increased, MR is gradually suppressed 
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FIG. 1. (a) Powder XRD patterns and structural refinement 
results. The data were shown by (+) , and the fit is given 
by the red solid line. The difference curve (the green solid 
line) is offset. The inset shows the single crystal XRD pat- 
tern indicating the c-axis orientation of crystal, (b) Tem- 
perature dependence of the in-plane resistivity pab{T) (open 
symbols) and c-axis resistivity pciT) (filled symbols) in B = 
T (squares) and B = 9 T (circles) magnetic field respectively. 



and is rather small above ~ 60 K. Resistivity along the 
c-axis (pc{T)) is nearly two orders of magnitude larger 
than pab{T) and exhibits a weak crossover at high tem- 
perature. In what follows we will discuss in-plane MR. 

Angular dependent MR 6*) at T ~ 2 K is shown 
in Fig. 2(a)-(b). The crystal was mounted on a rotating 
stage such that the tilt angle between sample surface 
(a6-plane) and the magnetic field can be continuously 
changed, with currents flowing in the a6-plane perpen- 
dicular to magnetic field (inset in Fig. 2(a)). The mag- 
netoresistancc of SrMnBi2 exhibits significant angular de- 
pendence (Fig. 2(a,b)). When B is parallel to the c-axis 
(6* = 0°), the MR is maximized and is linear above a char- 
acteristic field (~' 1 T). With increase in the tilt angle 6*, 
MR gradually decreases and becomes nearly negligible 
for B in the afo-plane {9 = 90°). 

The response of the carriers to the applied mag- 
netic field and the magnitude of MR is determined by 
the mobility in the plane perpendicular to the mag- 
netic fieldiii For nearly isotropic three-dimensional (3D) 
FS, there should be no significant angle-dependent MR 
(AMR). In (quasi-)2D systems, 2D states will only 
respond to perpendicular component of the magnetic 
field i3|cos(6')| and consequently longitudinal AMR and 
AMR-oseillation were observed in some quasi-2D conduc- 
tors (such as Sr2Ru04 and /3-(BEDT-TTF)2l3) and the 
surface states of TlsM^^^^LB. For example, the MR of 
the bulk state in topological insulator has only ^ 10% an- 
gular dependence while the angular dependence of MR in 
the surface state is about ten times larger Significant 




16 20 24 28 32 36 10 20 30 40 
B (T) Tilt angle 9 (Degree) 



FIG. 2. (a) In-plane resistivity p vs. the tilt angle 9 from 
0° to 360° at B = 3, 6 and 9 T and r = 2 K. The red line 
is the fitting curve using | cos(&)| (see text). Inset shows the 
configuration of the measurement, (b) In- plane Resistivity p 
vs. magnetic field B with different tilt angles 6 at 2 K. (c) 
MR SdH oscillations Ai? = Rxx ~ < R > as a function of field 
B below 35 T with tilt angles 6 from 0° to 35° at 2 K. The 
dashed lines indicate the SdH dips at Landau filling factor 
V — 5 and 6. The different curves are offset for clarification, 
(d) Position of dips with v — 5 and 6, as well as the position 
of peaks with v = 5, plotted against the tilt angle 9. The data 
are consistent with the 1/| cos(S)| dependence (red lines). The 
inset shows the Fourier transform of the SdH oscillation which 
gives a single frequency of _F = 138(7) T. 



AMR was also observed in some materials with highly 
anisotropic 3D FS such as Bi and Cu, but the period of 
AMR is determined by the shape of the Fermi surface and 
is very different from the one in 2D systems. In Bi, elec- 
trons exhibit a threefold valley degeneracy and in-plane 
mass anisotropy, so the AMR peaks each time when the 
magnetic field is oriented along the bisectrix axis and has 
a 60° periodi^ In Cu the AMR is more complex due to 
the complex FSs and peaks about every 25° 

The electronic structure calculations in SrMnBi2 show 
that the states near the Fermi energy Ep are dominated 
by the Bi states in the Bi square nets. Consequently the 
dominant FS should be quasi-2D. Angular dependent re- 
sistivity in i? = 9 T and T = 2 K shows wide maximum 
when the field is parallel to the c-axis {9 — 0°, 180°), and 
sharper minimum around 9 = 90°, 270° (Fig. 2(a)). The 
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whole curve of AMR in SrMnBi2 follows the function of 
|cos(0)| very well with a 180° period (red line in Fig. 
2(a)). Moreover, the larger pc than pab in Fig. 1(b) im- 
plies that the transfer integral and the coupling between 
layers along the c-axis is very small. All this implies that 
the FSs in SrMnBi2 should be highly anisotropic and that 
the mobility of carriers along is much smaller than the 
value in k^ky plane. 

Angular dependent MR quantum oscillations arc di- 
rectly related to the cross section of FS. In Fig. 2(c), 
the in-plane Ai? = R— < R > measured using same 
configuration as shown in the inset of Fig. 2(a) exhibits 
clear Shubnikov-de Hass (SdH) oscillations with tilt an- 
gles 6 from 0° to 35°. In metals, SdH oscillations cor- 
respond to successive emptying of LLs by the magnetic 
field and the LLs index n is related to the cross section of 
FS Sf by 27r(n+7) = 5i.^,iiii^>^ For a 2D FS (a cylin- 
der), the cross section has Sf{0) = 5*0/1 cos(0)| angular 
dependence and the LLs positions should be inversely 
proportional to |cos(6*)|j ^^'^^ The peak (dip) positions in 
SrMnBi2 rapidly shift toward higher field direction with 
increase in 9 (as indicated by the dashed lines in Fig. 
2(c)). In Fig. 2(d), the dip positions corresponding to 
LLs n = 5, 6 and the peak position with n = 5 were 
plotted against the tilt angle and can be described very 
well by 1/1 cos(6')| (the red lines in Fig. 2(d)). Similar be- 
havior was observed in the surface states of Tl^^^ and 
some other layered structures i^i^^ The Fourier trans- 
form of the SdH oscillation (inset of Fig. 2(d)) 
revealed that the oscillation component shov^rs a 
periodic behavior in 1/i? with a single frequency 
F = 138(7) T. The small value of frequency is con- 
sistent with previous value in Ref. [7] and demon- 
strates that the dominant FSs are very small since 
the Onsager relation is F = ($o/27r^)ylfc where $o 
is the flux quantum and Ak is the cross sectional 
area of FS^^ This clearly shows that the domi- 
nant two-dimensional FSs found in Fig. 2(b) are 
indeed the small FSs between F and M points, 
rather than the large FS at F point in the Bril- 
louin zone. Above SdH oscillation combined with the 
angular MR clearly suggests that the dominant FSs of 
SrMnBi2 are small quasi-2D cylinders along kz originat- 
ing from Bi square nets. In addition, there are still con- 
ventional parabolic bands with three dimensional char- 
acteristic close to the Fermi level^i causing the small de- 
viation from quasi-2D transport. 

Now we turn to the linear nonsaturated in-plane mag- 
netoresistancc in SrMnBi2 (Fig. 3(a)). The MR is linear 
over a wide field range up to 50 K. This behavior ex- 
tends to very low fields where the MR naturally reduces 
to a weak-field semiclassical quadratic dependence. The 
cross-over from the weak-field dependence to the high- 
field linear dependence can be best seen by considering 
the field derivative of the MR, dMR/dB (Fig. 3(b)). In 
the low field range [B <1 T at 2 K), dMR/dB is pro- 
portional to B (as shown by lines in low-field regions), 
indicating the semiclassical MR ^ A2B^. But above a 
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FIG. 3. (a) The magnetic field (B) dependence of the in- 
plane magnetoresistance (MR = {p{B) — p(0))/p(0)) at dif- 
ferent temperatures, (b) The field derivative of in-plane MR 
at different temperature respectively. The lines in high field 
regions were fitting results using MR = AiB + 0{B^) and the 
lines in low field regions using MR — A2B^. 




r(K) 



FIG. 4. (a) Temperature dependence of the critical field B* 
(black squares) and the high field MR linear coefficient Ai 
(blue circles) up to 50 K. The red solid line is the fitting results 
using B* — 2eBi,ii {EF + ksT)'^. (b) Temperature dependence 
of the effective MR mobility pmr extracted from the weak- 
field MR. 
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characteristic field B* , dMK/dB deviates from tfie semi- 
classical behavior and saturates to a much reduced slope 
(as shown by lines in the high- field region). This indi- 
cates that the MR for B > B* is dominated by a lin- 
ear field dependence plus a very small quadratic term 
(MR= AiB + 0{B^)). 

The linear MR deviates from the scmiclassical B^ de- 
pendence of MR in the low field region and a satu- 
rating MR in high fieldsi^i The unusual nonsaturating 
linear magnetoresistance has been reported in gapless 
semiconductor Ag2-g(Te/Se) ^''i^^ with the linear energy 
spectrum in the quantum limit J^i^ Recent first princi- 
ple calculations confirmed that these materials have a 
gapless Dirac-type surface state.— Linear magnetoresis- 
tance is also observed in topological insulatorai^*^ and 
BaFe2As2i^ with Dirac fcrmion states. Another possible 
origin of the large linear magnetoresistance is the mobil- 
ity fluctuations in a strongly inhomogeneous system)^ 
This does not apply in SrMnBi2 since our sample is sto- 
ichiometric crystal without doping/disorder. Below we 
show that the nonsaturating linear magnetoresistance 
and the deviation from the scmiclassical transport in 
SrMnBi2 is due to the linear energy dispersion. 

The application of a strong perpendicular external 
magnetic field (B) would lead to a complete quan- 
tization of the orbital motion of carriers with linear 
energy dispersion, resulting in quantized LLs £'„ = 
sgn{n)v F 2ehB\n\ where n = 0,±1,±2,--- is the LL 
index and is the Fermi velocity] ^^'^'^ Then the energy 
splitting betw een th e lowest and 1^* LLs is described by 
All ~ ±vp\/2ehBJ^^ In the quantum limit at specific 
temperature and field, All becomes larger than both the 
Fermi energy Ep and the thermal fiuctuations kgT at a, 
finite temperature. Consequently all carriers occupy the 
lowest Landau level and eventually the quantum trans- 
port with linear magnetoresistance shows up. The crit- 
ical field B* above which the quantum limit is satisfied 
at specific temperature T is B* = 2^^^{Ef + ksT)^ 
The splitting of LLs in conventional parabolic bands is 
All ~ Hence the evolution of All with field 

for parabolic bands is much slower than that for Dirac 
fermion states, and it is difficult to observe quantum limit 
behavior in the moderate field range. The temperature 
dependence of critical field B* in SrMnBi2 clearly devi- 
ates from the linear relationship and can be well fitted by 
B* = ^^{EF + kBT)"^, as shown in Fig. 4(a). The fit- 
ting gives the Fermi velocity vp ^ 5.13 x 10^ ms~^ and 
Ai ~ 4.97 mcV. This confirms the existence of Dirac 



fermion states in SrMnBi2. 

In a multiband system with both Dirac and conven- 
tional parabolic-band carriers (including electrons and 
holes) the magnetoresistance in the scmiclassical trans- 
port can be described as MR = ^^^^^^f^^p^^^ where 
ae,crh-, fJ-e, tJ-h are the effective electron and hole conduc- 
tivity and mobility in zero field respectively, when the 
Dirac carriers are dominant in transport j ^^i^^ Then the 
coefficient of the low-field B^ quadratic term, A2, is re- 
lated to the effective MR mobility y/A^ = ^^^^ (A*e + 
/i/i) ~ fJ-MR^ which is smaller than the average mobil- 
ity of carriers Have = and gives an estimate of 
the lower bound. Fig. 4(b) shows the dependence of 
fiMR on the temperature. At 2 K, the value of ^mr 
is about 3400 cm^/Vs. The large effective MR mobility 
also implies that Dirac fermions dominate the transport 
behavior. With increase in temperature, the value of 
/iA/fl and the coefficient of high- field linear term Ai (Fig. 
4(a)) decrease sharply. This is due to thermal fiuctuation 
smearing out the LL splitting. 

In summary, we demonstrate quantum transport of 
2D Dirac fermion states in bulk SrMnBi2 single crys- 
tals. The bands with linear energy dispersion lead to 
the large nonsaturated linear magnetoresistance since all 
Dirac fermions occupy the lowest Landau level in the 
quantum limit. The transverse magnetoresistance ex- 
hibits a crossover at a critical field B* from scmiclassi- 
cal weak-field B^ dependence to the high-field linear-field 
dependence. With increase in temperature, the critical 
field B* increases and the temperature dependence of B* 
satisfies quadratic behavior which is attributed to the 
Landau level splitting of the linear energy dispersion. 
The effective magnetoresistant mobility fiMR ~ 3400 
cm^/Vs comparable to values observed in graphene is 
observed. The angle-dependence of magnetoresistance 
shows a |cos(6')| dependence while the LL positions in 
SdH oscillations are inversely proportional to |cos(6')|, 
indicating the dominant quasi-2D Fermi surfaces. Our 
results show that the crystals with Bi square nets can 
host phenomena commonly observed so far in 2D struc- 
tures and materials like graphene. 
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